Analysis of a quenched lattice-QCD dressed-quark propagator 
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Quenched lattice-QCD data on the dressed-quark Schwinger function can be correlated with 
dressed-gluon data via a rainbow gap equation so long as that equation's kernel possesses enhance- 
ment at infrared momenta above that exhibited by the gluon alone. The required enhancement 
can be ascribed to a dressing of the quark-gluon vertex. The solutions of the rainbow gap equa- 
tion exhibit dynamical chiral symmetry breaking and are consistent with confinement. The gap 
equation and related, symmetry-preserving ladder Bethe-Salpeter equation yield estimates for chi- 
ral and physical pion observables that suggest these quantities are materially underestimated in the 
quenched theory: |(gg)| by a factor of two and by 30%. 

PACS numbers: 12.38.Lg, 12.38.Gc, 12. 38. Aw, 24.85. +p 



I. INTRODUCTION 

It is a longstanding prediction that the Schwinger func- 
tions which characterise the propagation of QCD's el- 
ementary excitations: gluons, ghosts and quarks, are 
strongly modified at infrared momentum scales, namely, 
spacelike momenta P < 2GeV^ Q,I1I3. Indeed, this 
property of asymptotically free theories was elucidated 
in Refs. 4, 5] and could be anticipated from studies of 
strong coupling QED Such momentum-dependent 
dressing is a fundamental feature of strong QCD that 
is observable in hadronic phenomena For example: 
it is the mechanism by which the current-quark mass 
evolves to assume the scale of a constituent-quark mass 
at infrared momenta, and thereby that dynamical chiral 
symmetry breaking (DCSB) is exhibited; and it may also 
provide an understanding of confinement, as we canvass 
in Sec. IIIIDl These are keystones of hadron physics |^. 

Numerical simulations of lattice-QCD provide direct 
access to QCD's Schwinger functions, and recent studies 
of the quenched theory yield dressed-gluon 0, 0, 
and -quark 0, 0| two-point functions ( "propagators" ) 
that are in semi-quantitative agreement with Dyson- 
Schwinger equation (DSE) calculations | il5> . il6i . .l7i.il8i . il9k 
123,1^ ■ However, these dressed-gluon and -quark propa- 
gators are not obviously consistent with each other in the 
following sense: use of the lattice dressed-gluon two-point 
function as the sole basis for the kernel of QCD's gap 
equation cannot yield the lattice dressed-quark propa- 
gator without a material infrared modification (enhance- 
ment) of the dressed-quark-gluon vertex [2^l23Ll2^ . For- 
tunately, just such behaviour can be understood to arise 
owin g to multiplicative renormalisability of the gap equa- 
tion j25l, (2^ and is observed in lattice estimates of this 
three-point function ^27, 28, 29]. 

Herein we elucidate these points using a concrete 
model for the gap equation's kernel. Naturally, since the 
ultraviolet behaviour of this kernel is fixed by perturba- 
tive QCD and hence model- independent, our study will 
focus on and expose aspects of the infrared behaviour of 



the Schwinger functions described above. 

Furthermore, as indicated at the outset, DCSB is en- 
coded in the chiral-limit behaviour of the dressed-quark 
propagator. However, contemporary lattice-QCD simu- 
lations are restricted to current-quark masses that are too 
large for unambiguous statements to be made about the 
magnitude of this effect. With a well-constrained model 
for the gap equation's kernel it is straightforward to cal- 
culate the dressed-quark propagator in the chiral limit. 
Hence our analysis will also provide an informed estimate 
of the chiral limit behaviour of the lattice results. 

Finally, we consider two additional questions; namely, 
how do Schwinger functions obtained in simulations of 
quenched lattice-QCD differ from those in full QCD, 
and can that difference be used to estimate the effect 
of quenching on physical observables? Our model for the 
gap equation's kernel provides a foundation from which 
we believe these problems can fruitfully be addressed. 

The article is organised as follows. In Sec^Jwe review 
the gap equation, the form of its solution and the nature 
of its kernel. Section Hill describes the construction of a 
model for the gap equation's kernel that correlates lattice 
results for the dressed-gluon and -quark two point func- 
tions, and explains aspects of the dressed-quark function 
obtained therewith. The kernel is exploited further in 
Sec. IIVI wherein it provides the basis for calculating in- 
formed estimates of observable quantities in quenched- 
QCD. Sectionals an epilogue. 



II. GAP EQUATION 



The renormalised dressed-quark propagator, S{p), is 
the solution of the DSE 



S-\p)^Z2{<:, A) Z7 • p -f Z4(C, A) m(C) + S'(p, A), (1) 
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wherein the dressed-quark self-energy is ^ 



E' [p, A) = Zi (C, A)^ g'D^, {p - q) ^J, S{q) K (g, p) 



Here D^^{k) is the renormalised dressed gluon propa- 
gator; T]j(q,p) is the renormalised dressed quark-gluon 
vertex; and Zi(^,A), Z2(C,A) and Zn{C^,K) are, re- 
spectively, Lagrangian renormalisation constants for the 
quark-gluon vertex, quark wave function and current- 
quark mass. 

In Eq. 121, :— d^q/ {2it)'^ denotes a translation- 
ally invariant ultraviolet regularisation of the momen- 
tum space integral, with regularisation mass-scale A. In 
practice a Pauli-Villars scheme is used, in which an ad- 
ditional factor of h?/{h? + [p — qY) is included with 
the gluon propagator. The resulting integral is finite 
VA < oo, and develops a logarithmic divergence when 
the regularisation is removed; i.e., A — > oo, which is 
the final stage of any calculation. We emphasise that 
it is only with a translationally invariant regularisation 
scheme that Ward-Takahashi identities can be preserved, 
something that is crucial to ensuring, e.g., axial- vector 
current conservation in the chiral limit. 

The general form of the dressed-quark propagator is 



S[p) 



z{p'-e 



■ p + M{p'^) 



(3) 



where the Lorentz scalar functions Z{p'^;('^) and M{p^) 
are respectively referred to as the quark wave function 
renormalisation and running quark mass (or dressed- 
quark mass-function). The renormalisation constants Z2 
and Zi are determined by solving the gap equation, Eq. 

subject to the renormalisation condition that at some 
large spacelike 



S-\p) 



i-/-p + m{C), 



(4) 



where to(C) is the renormalised current-quark mass. The 
renormalised dressed-quark propagator is independent of 
the regularisation mass-scale, A. It depends on the renor- 
malisation point, in a manner prescribed by the the- 
ory's dynamics. This dependence is expressed in the cal- 
culable ^-dependence of the wave function renormalisa- 
tion. The dressed-quark mass-function is independent 
of the regularisation mass-scale and of the renormalisa- 
tion point. Once the renormalisation scheme has been 
faithfully applied the regularisation mass-scale may be 
removed to infinity. 

Given the form of the dressed gluon propagator, 
Dp,i,{p — q), and dressed quark-gluon vertex, r*j(g,p), it is 



straightforward to determine the corresponding dressed 
quark propagator using well-established numerical meth- 
ods. A particularly important characteristic of the non- 
linear integral equation in Eq. ^ is that it yields a 



^2) nonzero solution for M{p'^) in the chiral limit 



t yields 

EE 113: 



Z4(C,A)m(C) =0, A»C: 



(5) 



if, and only if, there is sufficient infrared support in the 
integrand. This is dynamical chiral symmetry breaking, 
which we discuss further in Sec. IIV Al It is noteworthy 
that for finite C, and A 00, the left hand side (l.h.s.) of 
Eq. is identically zero, by definition, because the mass 
term in QCD's Lagrangian density is renormalisation- 
point-independent. The condition specified in Eq. jSJ, 
on the other hand, effects the result that at the (per- 
turbative) renormalisation point there is no mass-scale 
associated with explicit chiral symmetry breaking, which 
is the essence of the chiral limit. 

As will now be plain, the kernel of the gap equation, 
Eq. (Q) , is formed from the product of the dressed-gluon 
propagator and dressed-quark-gluon vertex. The equa- 
tion is therefore coupled to the DSEs satisfied by these 
functions. Those equations in turn involve other n-point 
functions and hence a tractable problem is only realised 
once a truncation scheme is specified. At least one non- 
perturbative, chiral symmetry preserving truncation ex- 
ists 1^ and the first term in that scheme is the 
renormalisation-group-improved rainbow gap equation, 
wherein the self-energy, Eq. (0), assumes the form 



J a 



(6) 



where Q = p — q and Di^^°{Q) is the Landau gauge free- 
gluon propagator. 

In Eq. ©, G{Q^) is an effective interaction, which ex- 
presses the combined effect of dressing both the gluon 
propagator and quark-gluon vertex consistent with the 
constraints imposed by, e.g., vector and axial- vector 
Ward-Takahashi identities. Asymptotic freedom entails 



(7) 



viz., the effective interaction is proportional to the strong 
running coupling in the ultraviolet. However, its form is 
unknown for ^ 2 GeV^. An explanation of many di- 
verse hadron phenomena has been obtained by modelling 
this behaviour Tj but our goal is different. Hereinafter 
we explore the ramifications of employing the dressed- 
gluon propagator inferred from numerical simulations of 
lattice-QCD in building the effective interaction. 



III. ANALYSING LATTICE DATA 



A. Lattice gluon propagator 

'l the Landau gauge dressed glu' 
that obey {7^,7^} = 25^^,. tor was computed using quenched lattice-QCD configu- 



[1] We use a Euclidean metric wherewith the scalar product of two 

fourvectorsisa ■ b = ^^^^ Oifci, andHermitian Dirac-7 matrices R-^f. ^ the Landau gauge dressed gluon propaga- 
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rations and the result was parametrised as: 



(fc2 + A2)l+" 



(8) 



with 



A 



9.i 



j+O.l 
'-0.9' 



a 



n 9+0.1+0.2 

^•^-0.2-0.3' 



Ag ^ 1.020 ±0.1 ±0.025 GeV, 

(9) 

where the first pair of errors are statistical and the sec- 
ond, when present, denote systematic errors associated 
with finite lattice spacing and volume. In the simulation 
the lattice spacing a — 1/[1.885 GeV]. The numerator in 
the second term of Eq. JHJ is 



L{k',Ag)^[ - In 



(10) 



with do = [39 - 9^ - 47V/]/[2(33 - 2iV/)], an expres- 
sion which ensures the parametrisation expresses the cor- 
rect one-loop behaviour at ultraviolet momenta. In the 
quenched, Landau-gauge study, Nf = , ^ = 0, so 



do = 13/22. 

B. Effective quark-gluon vertex 

We have noted that for Q'^ > 2 GeV^ in Eq. © 



(11) 



ln(gVA2 



(12) 



QCD^ 



where jm = 12/(33 — 2A^/) is the anomalous mass di- 
mension. To proceed we therefore write 



^giQ')^DiQ')T,iQ'), 

with D{Q^) given in Eq. © and 

1 [iln(r + QVA^)]'^° 



Zg [ln(T + 02/A2 



QCD 



)] 



(13) 



«(Q'), (14) 



where r = — 1 > 1 is an infrared cutoff. Equation 
(|13|l factorises the effective interaction into a contribu- 
tion from the lattice dressed-gluon propagator multiplied 
by a contribution from the vertex, which we shall subse- 
quently determine phenomenologically. 

We remark that the renormalisation-group-improved 
rainbow truncation retains only that single element of 
the dressed-quark-gluon vertex which is ultraviolet di- 
vergent at one-loop level and this explains the simple 
form of Eq. (|14() . Systematic analyses of corrections to 
the rainbow truncation show Fi to be the dominant am- 
plitude of the dressed vertex: the remaining amplitudes 
do not significantly affect observables [32] . In proceeding 



phenomenologically solely with Fi, we force f (Q^) to as- 
sume the role of the omitted amplitudes to the maximum 
extent possible. 

In Eq. (ini), so long as v{Q^) ~ 1 for > 2 GcV^, Eq. 
(|12|l is satisfied and consequently the rainbow gap equa- 
tion preserves the renormalisation group fiow of QCD at 
one-loop. We therefore consider a simple Ansatz with 
this property: 



vm = 



a^{m) + Q^/Al 
h^Q^Al 



where 

ay{m) 



ai 



1 + 02 [m{()/Ag\ + 03 [m{()/Ag 



(15) 



(16) 



and ai_2.3 and b are dimensionless parameters, which are 
fitted by requiring that the gap equation yield a solution 
for the dressed-quark propagator that agrees well point- 
wise with the results obtained in numerical simulations 
of quenched lattice-QCD [ll[l|. It is important to note 
that a good fit to lattice data is impossible unless (m) 
depends on the current-quark mass. While more compli- 
cated forms are clearly possible, the Ansatz of Eq. H16|) 
is adequate. 



C. Fit to lattice results 

Now, to be explicit, the parameters in Eqs. ((T3|l . 
were determined by the following procedure. The rain- 
bow gap equation; viz., Eq. |^ simplified via Eq. JHJ, 
was solved using the effective interaction specified by Eqs. 
(inil-llini), with D{Q'^) exactly as given in Eq. ijSJi. 

The ultraviolet behaviour of the mass function, M (p^), 
is determined by perturbative QCD and is therefore 
model independent. Hence the current-quark mass, 
m(C), was fixed by requiring agreement between the DSE 
and lattice results for M{p'^) on > IGeV^. We se- 
lected three lattice data sets from Ref. Q and, for con- 
sistency with Refs. 0,0], used a renormalisation point 
C = 19 GeV, which is well into the perturbative domain. 
This gave 



a Wlatticc 



m(C)(GeV) 



0.018 0.036 0.072 



0.030 0.055 0.110 



(17) 



The dimensionless parameters ai^2,3 and 6 were subse- 
quently determined in a simultaneous least-squares fit of 
DSE solutions for M{p'^) at these current-quark masses 
to all the lattice data. This necessarily required the gap 
equation to be solved repeatedly. Nevertheless, the fit re- 
quired only hours on a modern workstation, and yielded: 



(18) 



ai 


02 


0-3 


b 


1.5 


7.35 


63.0 


0.005 



These parameters completely determine the "best-fit 
effective-interaction" and hence our lattice-constrained 
model for the gap equation's kernel. 
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FIG. 1: Data, upper three sets: lattice results for M{p'^) in 
GeV at am values in Eq. (1171 : lower points (boxes): linear 
extrapolation of lattice results ,14] to am — 0. Solid curves: 
best-fit-interaction gap equation solutions for M{p^) obtained 
using the current-quark masses in Eq. ; dashed-curve: gap 
equation's solution in the chiral limit, Eq. ©. 

We emphasise that because the comparison is with 
simulations of quenched lattice-QCD, we used Nf = 
throughout and jHlliliilll 

Aqu-QCD = 0.234 GeV. (19) 

The strength of the running strong coupUng is underesti- 
mated in simulations of quenched lattice-QCD 0. (NB. 
Halving or doubling Aqu_QCD has no material quantita- 
tive impact on the results reported in Sec. IIVI nor does 
it qualitatively affect our conclusions.) 



1. Fidelity and quiddity of the procedure 

In Fig. n] we compare DSE solutions for A/(p^), ob- 
tained using the optimised effective interaction, with lat- 
tice results. In addition, we depict the DSE solution 
for M{p^) calculated in the chiral limit along with the 
linear extrapolation of the lattice data to am ~ 0, as 
described in Ref. It is apparent that the lattice- 

gluon and lattice-quark propagators can be correlated 
via the renormalisation-group-improved gap equation. 
That was achieved via v{Q^) in Eq. IjlSI) . and the re- 
quired form is depicted in Fig. [3 Plainly, consistency 
between the propagators via this gap equation requires 
an infrared enhancement of the vertex, as anticipated 
in Refs. l25l . Our inferred form is in semi- 

quantitative agreement with the result of recent, ex- 
ploratory lattice-QCD simulations of the dressed-quark- 
gluon vertex [27Ll2ll29t. 

Dynamical chiral symmetry breaking is another im- 
portant feature evident in Fig. ^ viz., the existence of a 
M(p^) ^ solution of the gap equation in the chiral limit. 
We deduce that DCSB is manifest in quenched-QCD and. 
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FIG. 2: Dimensionless vertex dressing factor: v{Q^), defined 
via Eqs. I|15^ . 116|l . CHJi obtained in the chiral limit {solid 
curve) and with the current-quark masses in Eq. (I17II . v{Q^) 
is finite at = and decreases with increasing m(("). 



in the following, quantify the magnitude of that effect. It 
should be observed that a linear extrapolation to am — 
of the lattice data obtained with nonzero current-quark 
masses overestimates the mass function calculated di- 
rectly as the solution of the gap equation. 

Figure 13 focuses on the lattice simulations for the 
intermediate value of the current-quark mass; namely, 
am = 0.036, and compares lattice output for both Af(p^) 
and the quark wave function renormalisation, Z(p^), with 
our results. We emphasise that the form of Z{j?) was not 
used in fitting f (Q^). Hence the pointwise agreement be- 
tween the gap equation's solution and the lattice result 
indicates that our simple expression for the effective in- 
teraction captures the dominant dynamical content and, 
in particular, that omitting the subdominant amplitudes 
in the dressed-quark-gluon vertex is not a serious flaw in 
this study. 



D. Spectral properties 

In a quantum field theory defined by a Euclidean mea- 
sure 'Sal the Osterwalder-Schrader axioms [s^ |33 are 
five conditions which any moment of this measure (n- 
point Schwinger function) must satisfy if it is to have 
an analytic continuation to Minkowski space and hence 
an association with observable quantities. One of these 
is "0S3" : the axiom of reflection positivity, which is vi- 
olated if the Schwinger function's Fourier transform to 
configuration space is not positive definite. The space of 
observable asymptotic states is spanned by eigenvectors 
of the theory's infrared Hamiltonian and no Schwinger 
function that breaches 0S3 has a correspondent in this 
space. Consequently, the violation of 0S3 is a sufficient 
condition for confinement. 

This connection has long been of interest [3^ , and 
is discussed at length in Refs. pol llll , and reviewed 
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FIG. 3: Data, quenched lattice-QCD results for M{p^) and 
Z(p^) obtained with am = 0.036 0|; dashed curve, Z{p^), 
and solid curve, M{p^) calculated from the gap equation with 
our optimised effective interaction and m((") — 55MeV. (NB. 
Z{p^) is dimensionless and M{p'^) is measured in GeV.) 
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FIG. 4; |Ag(r)| obtained from: the chiral limit gap equation 
solution calculated using our lattice-constrained kernel, solid 
curve; Eq. with ct = 0.13 GeV, 9 = 7r/2.46, dotted curve; 
the model of Ref. dashed curve. 



in Sec. 6.2 of Ref. Q, Sec. 2.2 of Ref. 1.] and Sec. 2.4 of 
Ref. l3|. It suggests and admits a practical test "2^ that 
has been exploited in Refs. 20, ,22, 23, ,43, 44, 45, 4fi,il 
|4^ figf and which for the quark 2-point function is based 
on the behaviour of 



As(r) 



(fx 



(27r)4 
1 1'°° 

- / de cos{eT)as{e^) 



(20) 
(21) 



where as is the Dirac-scalar projection of the dressed- 
quark propagator. For a noninteracting fermion with 
mass /i, 



^frcc 



(T) 



1 

- / de cos(e T) 



(22) 



The r.h.s. is positive definite. It is also plainly related 
via analytic continuation (T it) to the free-particle 
solution of the Minkowski space Dirac equation. The 
existence of an associated asymptotic state is indubitable. 
If instead one encountered a theory in which [4lll43 . l49| 



1 



(23) 



+ — tp^ p^ + fi^ + tp^ 
a function with poles at p^ + exp(±i6') = 0, where 

cr* = // + p^, tan 9^ pV/^^ (24) 

then 



As(T) = ^e--^-^icosfaTsin| 



(25) 



This Fourier transform has infinitely many, regularly 
spaced zeros and hence 0S3 is violated. Thus the fermion 



described by this Schwinger function has no correspon- 
dent in the space of observable asymptotic states. 

It is readily apparent that Eq. (|23|l evolves to a free 
particle propagator when p 0. This limit is expressed 
in Eq. via 9 ^ 0, a ^ p, wherewith Eq. ||22Jl is 
recovered; a result that is tied to the feature that the 
first zero of As(T) in Eq. H25(l occurs at 



Zl 



2(7 



(26) 



and hence zi — * cxo for p ^ |43l |44| . 

We calculated As{T) using the gap equation solutions 
discussed above and the form of |As(T)| obtained with 
the chiral limit solution is depicted in Fig. 01 The vi- 
olation of 0S3 is manifest in the appearance of cusps: 
ln|A5(r)| is negative and infinite in magnitude at zeros 
of As{T). The differences evident in a comparison with 
the result obtained from Eq. (|25|l indicate that the sin- 
gularity structure of the dressed-quark 2-point Schwinger 
function obtained from the lattice-constrained kernel is 
more complicated than just a single pair of complex con- 
jugate poles. However, the similarities suggest that this 
picture may serve well as an idealisation 49] . (NB. Qual- 
itatively identical results are obtained using (Ty(p^), the 
Dirac-vector projection of the dressed-quark propagator, 
instead of as{p^)) 

Figure 21 also portrays the result obtained with the ef- 
fective interaction proposed in Ref. jlaand used effica- 
ciously in studies of meson properties Significantly, 
the first zero appears at a smaller value of T in this case. 
Using the model of Eq. as a guide, that shift indi- 
cates a larger value of a. This sits well with the fact that 
the mass-scale generated dynamically by the interaction 
of Ref. ^3 is larger than that produced by the inter- 
action used herein, which is only required to correlate 



quenched lattice data for the gluon and quark two-point 
functions. 

From the results and analysis reported in this sub- 
section, we deduce that light-quarks do not appear in 
the space of observable asymptotic states associated with 
quenched-QCD; an outcome anticipated in Ref. p^ . 



IV. CHIRAL AND PHYSICAL PION 
OBSERVABLES 

A. Chiral limit 

The scale of DCSB is measured by the value of the 
renormalisation-point dependent vacuum quark conden- 
sate, which is obtained directly from the chiral limit 
dressed-quark propagator yO] : 

.A 

- {qq)°. = lim Z^{C,A) tr / So{q), (27) 



where "tr" denotes a trace only over Dirac indices and 
the index "0" labels a quantity calculated in the chiral 
limit, Eq. (|3J). In Eq. (|27|l the gauge parameter depen- 
dence of the renormalisation constant Z4 is precisely that 
required to ensure the vacuum quark condensate is gauge 
independent. This constant is fixed by the renormalisa- 
tion condition, Eq. (0J, which entails 



Z4(C,A) = - 



1 1 



m(C) 4 



-trS'(C,A) 



(28) 



wherein the right hand side is well-defined in the chi- 
ral limit Tt^. (It may also be determined by study- 
ing the fully amputated pseudoscalar-quark-antiquark 3- 
point function |3Clj|.) A straightforward calculation using 
our chiral limit result; i.e, the propagator corresponding 
to the dashed-curve in Fig. ^ yields 



C=19 GcV 



(0.22 GeV)' 



(29) 



To evolve the condensate to a "typical hadronic scale," 
e.g., C = 1 GeV, one may use [s^ 

(30) 

where Z„i is the gauge invariant mass renormalisation 
constant. Contemporary phenomenological approaches 
employ the one-loop expression for Z^ and following this 
expedient we obtain, practically as a matter of definition. 



mil GcV 



/ ln[l/Aqu^QCD] Y 

^ln[19/Aqu_QCD]y 

(0.19GeV)^ 



(-( 



99)?9Gov) 



(31) 



which may be compared with a best-fit phenomenological 
value j51,]: (0.24 ± 0.01 GeV)^. It is notable that DSE 
models which efficaciously describe light-meson physics; 
e.g., Refs. mill, give m^lccv = -(0.24 GeV)^. 
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FIG. 5: M{pf^ = 0.38 GeV^), in GeV, as a function of the 
current-quark mass. Solid curve, our result; circles, lattice 
data for am in Eq. (amvalues) fl4| : dashed-line, linear fit to 
the lattice data, Eq. 1321 . 



Our gap equation assisted estimate therefore indicates 
that the chiral condensate in quenched-QCD is a factor 
of two smaller than that which is obtained from analyses 
of strong interaction observables. These results are in 
quantitative agreement with Ref. 26]. 

A fit to the linear extrapolation of the lattice data; 
viz., to the boxes in Fig. gives a significantly larger 
value [ij: ~{qq)iGcV = (0.270 ± 0.027 GeV)^. How- 
ever, the error is purely statistical. The systematic error, 
to which the linear extrapolation must contribute, was 
not estimated. That may be important given the dis- 
crepancy, conspicuous in Fig.^ between our direct eval- 
uation of the chiral limit mass function and the linear 
extrapolation of the lattice data to am — 0: the linear 
extrapolation lies well above the result of our chiral limit 
calculation. 

This point may be illustrated further. In Fig. |S1 we 
plot M(p^ = Pjr), where piR = 0.62 GeV, as a function 
of the current-quark mass (solid curve). This value of 
the argument was chosen because it is the smallest 
for which there are two lattice results for M(p^) at each 
current-quark mass in Eq. H17() . Those results are also 
plotted in the figure. It is evident that on the domain of 
current-quark masses directly accessible in lattice simu- 
lations, the lattice and DSE results lie on the same linear 
trajectory, of which 



0.38 GeV^ m(C)) = 0.18 + 1.42 m(C) (32) 



provides an adequate interpolation. However, as appar- 
ent in the figure, this fit on am G [0.018, 0.072] provides 
a poor extrapolation to m{() = 0, giving a result 40% 
too large. 

In Fig. El we repeat this procedure, focusing solely on 
our value of M{p'^ = 0) because directly calculated lat- 
tice data are unavailable at this extreme infrared point 
and published estimates obtained by extrapolating func- 
tions fitted to the lattice-p^-dependence are inconsistent 
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0.125 



FIG. 6: Solid curve, calculated M{p^ = 0), in GeV, as a func- 
tion of the current-quark mass m{(). The circles mark the 
current-quark masses in Eq. Ill7t . Dashed-line, linear inter- 
polation of our result for M(p'^ — 0) on this mass domain. 



[l3|. The pattern observed in Fig. [SI is again visible. 
On the domain of current-quark masses for which lat- 
tice data are available, the mass-dependence of MijP) is 
well-approximated by a straight line; namely. 



M(p2 = 0, m(C)) = 0.37 + 0.68 m(C) ., 



(33) 



but the value of M(0) determined via extrapolation to 
to(C) = is 14% too large. 

In all cases our calculated result possesses significant 
curvature. At fixed p^, M(p^,m(C)) is a monotonically 
increasing function of m(^) but, while it is concave-down 
for m(C) ^ 0.1 GeV, it inflects thereafter to become 
concave-up. In addition, at fixed to(C), M{p^) is a mono- 
tonically decreasing function of p^. It follows that a linear 
extrapolation determined by data on am £ [0.018,0.072] 
will necessarily overestimate M{p'^) at all positive values 
of p^. The figures illustrate that the error owing to ex- 
trapolation increases with increasing p^ and hence a sig- 
nificant overestimate can be anticipated on the domain 
in which the condensate was inferred from lattice data. 

It is straightforward to understand the behaviour evi- 
dent in Figs, m qualitatively. The existence of DCSB 
means that in the neighbourhood of the chiral limit a 
mass-scale other than the current-quark mass determines 
the magnitude of M(p^). As the current-quark mass in- 
creases from zero, its magnitude will come to affect that 
of the mass function. The gap equation is a nonlinear 
integral equation and hence this evolution of the mass- 
dependence of M{p^) will in general be nonlinear. Only 
at very large values of the current-quark mass will this 
scale dominate the behaviour of the mass function, as 
seen in studies of heavy-quark systems j52| , and the evo- 
lution become linear. 



B. Pion properties 

The renormalised homogeneous Bethe-Salpeter equa- 
tion (BSE) for the isovector-pseudoscalar channel; i.e., 
the pion, is 



[nik;P) 



[xi{q;P)]srKi:iq,k;P), (34) 



where k is the relative momentum of the quark-antiquark 
pair, P is their total momentum; and 



xi(g;P)-5(<z+)ri(g;P)5(g_) 



(35) 



with r^^{k;P) the pion's Bethe-Salpeter amplitude, 
which has the general form 



ri(fc;P)=T^75 



iEAk;P)+i-PFAk;P) 



j-kk-PG^{k;P)+ a^, k,,P, H,{k- P) 



(36) 



In Eq. H34|) . K{q,k;P) is the fully-amputated quark- 
antiquark scattering kernel, and the axial-vector Ward- 
Takahashi identity requires that this kernel and that of 
the gap equation must be intimately related. The conse- 
quences of this are elucidated in Refs. j^H^I and in the 
present case they entail 



Ki:iq,k;P) = -g{ik-qr) 



"A" 






Y7. 


ts 





(37) 



which provides the renormalisation-group-improved lad- 
der-truncation of the BSE. The efficacy of combining 
the renormalisation-group-improved rainbow-DSE and 
ladder-BSE truncations is exhibited in Ref. 0- In par- 
ticular, it guarantees that in the chiral limit the pion is 
both a Goldstone mode and a bound state of a strongly 
dressed-quark and -antiquark, and ensures consistency 
with chiral low energy theorems [U . 

All the elements involved in building the kernel of the 
pion's BSE were determined in the last section and hence 
one can solve for the pion's mass and Bethe-Salpeter am- 
plitude. To complete this exercise practically we consider 



[i-e(p')] 



ri(fc;P) 



[xUq;P)]srKi:{q,k;P). 

(38) 

This equation has a solution for arbitrary P^ and solv- 
ing it one obtains a trajectory i{P'^) whose first zero 
coincides with the bound state's mass, at which point 
Ti^{k; P) is the true bound state amplitude. In general, 
solving for £{P'^) in the physical domain: P^ < 0, re- 
quires that the integrand be evaluated at complex values 
of its argument. However, as <C M^(0); i.e., the 
magnitude of the zero is much smaller than the charac- 
teristic dynamically-generated scale in this problem, we 
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Calc. (quenched) 


Experiment 




0.1385 


0.1385 


fn 


0.066 


0.0924 




0.063 





TABLE L Pion-related observables calculated using our 
lattice-constrained eflective interaction. m{(^ — 19 GeV) = 
3.3 MeV was chosen to give = 0.1395 GeV. The index "0" 
indicates a quantity obtained in the chiral limit. 

avoid complex arguments by adopting the simple expe- 
dient of calculating £{P'^ > 0) and extrapolating to lo- 
cate its timelike zero. The Bethe-Salpeter amplitude is 
identified with the = solution. Naturally, this ex- 
pedient yields the exact solution in the chiral limit; and 
in cases where a comparison with the exact solution has 
been made for realistic, nonzero light-quark masses, the 
error is negligible 55] , as we shall subsequently illustrate. 

Once its mass and bound state amplitude are known, 
it is strai ghtf orward to calculate the pion's leptonic decay 
constant jSfl: 

U S^^ 2P, = Z2 tr /^T'757MXi('Z; P) ■ (39) 

Jq 

In this expression, the factor of Z2 is crucial: it en- 
sures the result is gauge invariant, and cutoff and 
renormalisation-point independent. (The Bethe-Salpeter 
amplitude is normalised canonically 56J.) 

Table|l]lists values of pion observables calculated using 
the effective interaction obtained in Sec. IIII Cl bv fitting 
the quenched-QCD lattice data. (The nonzero current- 
quark mass in the table corresponds to a one-loop evolved 
value of m(lGeV) = 5.0 MeV.) To aid with the consid- 
eration of these results we note that unquenched chiral- 
limit DSE calcvilations that accurately describe hadron 
observables give ,1^J /° — 0.090 GeV. We infer from these 
results that the pion decay constant in quenched lattice- 
QCD is underestimated by 30%. Quantitatively equiva- 
lent results were found in Ref. . 

The rainbow-ladder truncation of the gap and Bethe- 
Salpeter equations is chiral symmetry preserving: with- 
out fine-tuning it properly represents the consequences of 
chiral symmetry and its dynamical breaking. The trun- 
cation expresses the model-independent mass formula for 
flavour-nonsinglct pseudoscalar mesons [soj a corollary of 
which, at small current-quark masses, is the Gell-Mann- 
Oakes-Renner relation: 

if^)'ml = -2m(C) + Oim'{C)) . (40) 

Inserting our calculated values on the l.h.s. and r.h.s. of 
Eq. igni), we have 

(0.093 GeV)^ cf. (0.091 GeV)^ (41) 

viz, the same accuracy seen in exemplary coupled DSE- 
BSE calculations; e.g., Ref. This establishes the 

fidelity of our expedient for solving the BSE. 



V. SUMMARY 

We studied quenched-QCD using a rainbow-ladder 
truncation of the Dyson-Schwinger equations (DSEs) and 
demonstrated that existing results from lattice simula- 
tions of qucnchcd-QCD for the dressed-gluon and -quark 
Schwinger functions can be correlated via a gap equation 
that employs a renormalisation-group-improved model 
interaction. As usual, the ultraviolet behaviour of this 
effective interaction is fully determined by perturbative 
QCD. For the infrared behaviour we employed an Ansatz 
whose parameters were fixed in a least squares fit of the 
gap equation's solutions to lattice data on the dressed- 
quark mass function, M(p^), at available current-quark 
masses. With our best-fit parameters the mass functions 
calculated from the gap equation were indistinguishable 
from the lattice results. The gap equation simultane- 
ously yields the dressed-quark renormalisation function, 
Z{p^), and, without tuning, our results agreed with those 
obtained in the lattice simulations. 

To correlate the lattice's dressed-gluon and -quark 
Schwinger functions it was necessary for the gap equa- 
tion's kernel to exhibit infrared enhancement over and 
above that observed in the gluon function alone. We at- 
tributed that to an infrared enhancement of the dressed- 
quark-gluon vertex. The magnitude of the vertex mod- 
ification necessary to achieve the correlation is semi- 
quantitatively consistent with that observed in quenched 
lattice-QCD estimates of this three-point function. 

With a well-defined effective interaction, the gap equa- 
tion provides a solution for the dressed-quark Schwinger 
function at arbitrarily small current-quark masses and, 
in particular, in the chiral limit: no extrapolation is in- 
volved. A kernel that accurately describes dressed-quark 
lattice data at small current-quark masses may therefore 
be used as a tool with which to estimate the chiral limit 
behaviour of the lattice Schwinger function. Our view is 
that this method is a more reliable predictor than a lin- 
ear extrapolation of lattice data to the chiral limit. Even 
failing to accept this perspective, the material difference 
between results obtained via the lattice-constrained gap 
equation and those found by linear extrapolation of the 
lattice data must be cause for concern in employing the 
latter. 

In addition, from a well-defined gap equation it is 
straightforward to construct symmetry-preserving Bethe- 
Salpeter equations whose bound state solutions describe 
mesons. We illustrated this via the pion, and calculated 
its mass and decay constants in our DSE model of the 
quenched theory. 

Finally, assuming that existing lattice-QCD data are 
not afflicted by large systematic errors associated with fi- 
nite volume or lattice spacing, we infer from our analysis 
that quenched QCD exhibits dynamical chiral symme- 
try breaking and dressed-quark two-point functions that 
violate reflection positivity but that chiral and physical 
pion observables are significantly smaller in the quenched 
theory than in full QCD. 
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